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ABSTRACT 

We report the detection of a double planetary system orbiting around the evolved 
intermediate-mass star HD 4732 from precise Doppler measurements at Okayama Astrophysi- 
cal Observatory (OAO) and Anglo- Australian Observatory (AAO). The star is a KO subgiant 
with a mass of 1.7 Mq and solar metallicity. The planetary system is composed of two giant 
planets with minimum mass of m sin i = 2.4 Af j, orbital period of 360.2 d and 2732 d, and eccen- 
tricity of 0.13 and 0.23, respectively. Based on dynamical stability analysis for the system, we 
set the upper limit on the mass of the planets to be about 28 Mj (i > 5°) in the case of coplanar 
prograde configuration. 

Subject headings: stars: individual: HD 4732 — planetary systems — techniques: radial velocities 



1. Introduction 

Precise radial-velocity measurements of stars 
have revealed more than 500 extrasolar plan- 
ets in the last 20 years. Thanks to the in- 
crease in velocity precision and the duration 
of the planet search programs, not only single 
planets but also many multiple-planet systems 
have been discovered. Around solar-type stars, 
including those systems with long-term radial- 



^ Department of Earth and Planetary Sciences, Tokyo 
Institute of Technology, 2-12-1 Ookayama, Meguro-ku, 
Tokyo 152-8551, Japan 

^Department of Astrophysics, School of Physics, Uni- 
versity of NSW, 2052, Australia 

^Okayama Astrophysical Observatory, National Astro- 
nomical Observatory of Japan, Kamogata, Okayama 719- 
0232, Japan 

*The Graduate University for Advanced Studies, 
Shonan Village, Hayama, Kanagawa 240-0193, Japan 

^National Astronomical Observatory of Japan, 2-21-1 
Osawa, Mitaka, Tokyo 181-8588, Japan 

® Hiroshima Astrophysical Science Center, Hiroshima 
University, Higashi-Hiroshima, Hiroshima 739-8526, Japan 

'^Nishi-Harima Astronomical Observatory, Center for 
Astronomy, University of Hyogo, 407-2, Nishigaichi, Sayo, 
Hyogo 679-5313, Japan 



velocity trends, ~30-50 % of stars with giant plan- 
ets a r e multiple-pla net systems (IFischer et al. 



ets a r e muitipie-pianet sy s tems 
I2OOII : IWright et al."1 120071 l2009h 



and systems 
comprising less massive planets such as super- 
Earths aiidJ^JegtunB-mas^ are more abun- 



Wittenmver et al.ll2011a ) 



dant dMavor et al. I l201lt [Howard et all l2012t 



The multiple-planet systems give us deep in- 
sights into the for mation and evolution of plane- 
tary systems fe.g. iFord 1120061 ). For example, evi- 



dence of planet-planet scattering events could be 
preserved in the current orbi tal configurations of 
planets (e.g. lFord et al.ll2005l ). and the prevalence 
of orbital mean-motion resonances constitutes 
strong support for differential co nvergent orbital 



migration having occu rred (e.g. iLee and Peale 



2002l : lKlev et al. 1120041 ). Their statistical proper- 



ties may also provide us hints about the mecha- 
nisms of these processes; the rather uniform semi- 
major axis distribution without the pileup of hot 
Jupiters and the jump at 1 AU that are seen in 
single planets, and systematically l ower eccentric- 
ities compared to single planets (jWright et al. 
,2009,1 . 
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Studying the long-term orbital stability of 
multiple-planet systems is also important. Since 
the best fit orbital parameters derived from radial- 
velocity data do not necessarily guarantee that 
the resulting orbits are stable over long periods 
of time, dynamical stability analysis can check 
and further constrain the system parameters, 
such that the so lutions co nsistent w ith the data 
are stable (e.g. Lovis et al . 20 oj; Pepe et al. 



20071: iNiedzielski et al ll2009aHWright et al. 112009 : 



Wittcn mver et al.ll2012bl ). 



Johnson et al.l l2011bl: iRobertson et al.l |20~ : 



Here we report the detection of a double giant- 
planet system around the evolved intermediate- 
mass star HD 4732 (KO IV, M = 1.7 Mq), one 
of the 300 tar gets of the Okaya ma Planet Search 
Program (e.g. Isato et"al] l2012l ). This discovery 
was facilitated by joint precise radial- velocity ob- 
servations made from the Anglo- Australian Obser- 
vatory. Discoveries of planets and multiple-planet 
systems around such intermediate-mass stars have 
been rapidly growing in number during the past 
several years. Although the number of known 
multiple-planet systems around intermediate- 
mass stars is still small, it is interesting that 
many of them appear to be in mean-motion 
resonances such as 24 Sex (2:1; 



2011bh . HD 102272 (4:l:lNiedzie: 
V Op h (6:1; lOuirren bach et alT 
20m and BD20+2457 (3:2; 



Johnson et al 



ski et al. 2009a) 



miTsato et al 
Niedzielski et al 



2009bl ). Although detailed dynamical studies are 
required to determine w hether these sy stems are 
truly in resonance (e.g. [ Robertson et a l. 2012^ 



Wittenmver. et al. 2012cl ). the existence of such 



systems suggests that the planetary orbital migra- 
tion process for intermediate-mass stars is similar 
to that for solar-type ones. 

This paper is organized as follows. The stel- 
lar properties are presented in section [2] and the 
observations are described in section [H Orbital 
solution are presented in section U and results of 
line shape analysis are described in section [5j Dy- 
namical stability analysis is presented in section |6] 
and section [7] is devoted to summary. 

2. Stellar Properties 

HD 4732 (HR 2 28, HIP 38 34) is listed in the 
Hipparcos catalog (lESAlll997h as a KO giant with 
the apparent F-band magnitude F = 5.90 and the 



parallax tt = 17.70 ±0.99 mas, giving the distance 
of 56.5±3.2 pc and the absolute visual magni- 
tude My = -1-2.04 taking into account the correc- 
tion of interstellar e xtinction Ay =0.10 based o n 
Arenou et all (|l992l )'s table (|Takeda et al.ll2008l ). 



Atmospheric parameters (effective temperature 
Tcff, surface gravity log 5, micro-turbulent veloc- 
ity Vt^ and Fe abundance [Fe/H]) of all the tar- 
gets for Okayama Plan e t Sear ch Program were de- 
rived by iTakeda et al. I (I2OO8I) based on the spec- 
troscopic approach using the equivalent widths of 
well-behaved Fe I and Fe II lines of iodine- free stel- 
lar spectra. Details of the pr ocedure and re s ultant 
param eters are presente d in iTakeda et al.l ( 2002 ) 
and ,Takeda et all (|2008l ). For HD 4732, they ob- 
tained Teff = 4959 K, log5 (cgs) = 3.16, vt = 1.12 
km s~^, and [Fe/H]= -1-0.01 as well as projected 
rotational velocity of usinirot = 1-45 km s~^. 

With use of these atmospheric parameters to- 
get her with My, a bolometric correction based on 
the Kurucd (1993)'s theoretical calculation, and 
theore ti cal evolutiona r y trac ks of lLeieune fc Schaerer 



(120011 ). 

minosity 

L = 



Takeda et al 



( 2008f ) also determined lu- 
L and mass M for HD 4732 to be 
15.5 Lq and M = 1.74 Mq. The stellar 
radius R was estimated to be i? = 5.4 i?0 using 
the Stefan-Boltzmann relationship and the mea- 
sured L and Toff. Although the star is listed in 
the Hipparcos catalog as a giant star, the star is 
better regarded as a subgiant rather than giant 
based on the surface gravity and its position on 
the HR diagram (Fig [T]). The stellar properties 
are summarized in Table [1] 

Hipparcos observations revealed photometric 
stability for the star down to cthip = 0.006 mag. 
Furthermore, the star shows no significant emis- 
sion in the core of Ca II HK lines as shown in 
Figure [21 which suggests that the star is chromo- 
spherically ina c tive. For subgiants like HD 4732, 
IJohnson et aD (|2010h reported a typical radial- 
velocity jitter, which arises from a number of phe- 
nomena intrinsic to the star such as granulation, 
oscillation, and activity, of about 5 m s^^. The 
value is comparable to that of 6 m s~^ estimated 
for HD 4732 in section H 
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Observation 



3.1. OAO Observations 



Observations of HD 4732 at OAO were made 
with tire 1.88 nr telescope and the High Dis- 
persi on Echelle Spectrograph (HIDES; llzumiural 
19991 ) from August 2004 to January 2012. A 
slit width of the spectrograph was set to 200 
lim (0.76") corresponding to a spectral resolu- 
tion {R = A/AA) of 67000 by about 3.3 pix- 
els sampling. For precise radial-velocity mea- 
sure ments, we used an iodine absorption cell (I2 
cell; iKambe et al.l 120021 ) , which provides a fidu- 
cial wavelength reference in a wavelength range of 
5000-5800A. We have obtained 48 data points 
of HD 4732 with signal-to- noise ratio S/N=60- 
250 pix^-'^ by an exposure time of 900-1800 sec 
depending on the weather condition. The reduc- 
tion of echelle data (i.e. bias subtraction, fiat- 
fielding, scattered-light subtraction, and spectrum 
extraction) was performed using the IRAP0 soft- 
ware package in the standard way. 

For radial-velocity analysis, we modeled I2- 
superposed s tellar spectr a (star+ 12) by the method 
detailed in Sato et al. 1 (I2OO2I) and ISato et al 



( 2012 ) , whi ch is based on t he me thod bv lButler et al 



19961 ) and IValenti et al.l (|1995D . In the method. 



we model a star-|-l2 spectrum as a product of a 
high resolution I2 and a stellar template spectrum 
convolved with a modeled point spread function 
(PSF) of the spectrograph. We obtain the stellar 
spectrum by deconvolving a pure stellar spectrum 
with the spectrograph PSF estimated from an I2- 
superposed B-star spectrum. We have achieved a 
long-term Doppler precision of about 4-5 m s~^ 
over a time span of 9 years. Measurement error 
was derived from an ensemble of velocities from 
each of ~300 spectral chunks (each ^3A long) in 
every exposure. We listed the derived radial ve- 
locities for OAO data in Table [5] together with the 
estimated uncertainties. 

3.2. AAT Observations 

As HD 4732 is near the southern limit for OAO 
{6 ^ —25°), it is desirable to make observations 



^IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation, USA. 



from a more southerly site. To improve the phase 
coverage for HD4732, in 2010 September we be- 
gan ob serving it with the UCLES echelle spectro- 
graph ([Diego et al. I I1991I) at the 3.9-metre Anglo- 
Australian Telescope (AAT). UCLES achieves a 
resolution of 45,000 with a 1" slit. An iodine ab- 
sorption cell provides wavelength calibration from 
5000 to 6200 A. The spectrograph PSF and wave- 
length calibration are derived from the iodine ab- 
sorption lines em bedded on every pixel of the spec 



truin by the cell (jValenti et al.lll995t iButler et al 
19961 ). The result is a precision Doppler veloc- 



ity estimate for each epoch, along with an inter- 
nal uncertainty estimate, which includes the ef- 
fects of photon-counting uncertainties, residual er- 
rors in the spectrograph PSF model, and variation 
in the underlying spectrum between the iodine- 
free template, and epoch spectra observed through 
the iodine cell. The photon-weighted mid-time of 
each exposure is determined by an exposure meter. 
This technique has been successfully used at the 
AAT by the Anglo- Australian Planet Search (e.g. 
Tinney et al. 2001, O'Toole et al. 2009, Jones et 
al. 2010, Tinne y et al. 2011a) and the Pa n-Pacific 
Planet Search ( Wittenmver et al. 2011b[ ). All ve- 
locities are measured relative to the zero-point de- 
fined by the template observation. AAT/UCLES 
precis ion velocities are obtained using the Austral 
code (|Endl et al.ll2000( ). 

We have obtained 19 AAT observations of 
HD4732, and an iodine- free template observation 
was obtained on 2010 Oct 25. Exposure times 
ranged from 400 to 1200 s, with a resulting S/N 
of --100-200 per pixel each epoch. The AAT data 
span a total of 636 days, and have a mean internal 
velocity uncertainty of 5.0 ms~^. The data are 
given in Table |3l 

4. Orbit Fitting and Planetary Parame- 
ters 

For candidate multiple-planet systems, and for 
planet candidates with orbital periods near one 
year, it is critical to obtain the most complete 
possible phase coverage, and it is ideal to inde- 
pendently confirm the signal(s) from independent 
observatories. A recent example is HD 38283b, 
a 0.3 4MjuD planet with a period P — 363.2±1.6 
days ( Tinnev et al. 2011bl) . That planet required 
12 years of observations to confirm, and a ro- 
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bust detection was further aided by the fact that 
HD 38283 is circumpolar and so could be observed 
year-round (although at high airmasses below the 
pole) . Another example is HD 159868c, which has 
a peri od P = 352.3±1.3 days (jWittenmyer et al 



2012bl) . also worryingly close to one year. Again, 
more than 9 years of AAT data were needed for a 
secure detection, and the 352-day period was con- 
firme d independently with Keck data over a 4- year 
span (jWittenrnver et al. I l2012bl) . 

Five years of radial-velocity observations of 
HD 4732 at OAO revealed evidence of a signal 
with a period near one year (^338 days). Since 
this star is near the southern limit for OAO, it 
is only observable for 5 months of the year, re- 
sulting in persistent phase gaps which make the 
confirmation of such a period extremely difficult. 
AAT observations in 2010-2011 filled in the crit- 
ical phase gaps and confirmed the orbital period 
suggested by the OAO data. Continued observa- 
tions from both telescopes in 2011-2012 have also 
revealed a second velocity signal, first manifesting 
as a residual trend, and subsequently as a second, 
long periodicity. 

To fit the two Keplerian signals evident in the 
HD 4732 data, we first employed a genetic algo- 
rithm. This approach has proven useful in previ- 
ous work where it was necessary to fit highly un- 
certain orbits with lon g periods near the total d u- 
ration of observatio ns ( Wittenmver et al. 2012al lbl: 
Horner et al.ll2012l) . We ran the genetic algorithm 
for 10,000 iterations, each consisting of 1000 indi- 
vidual trial fits. The best-fit set of parameters 
is thus the result of ^ 10^ trial fits. The pa- 
rameters of the best 2-planet solution obtained 
by the genetic algorithm were the n used as ini 



tial inputs for the GaussFit code (jJeffervs et al 



119871 ). a nonlinear least-squares fitting routine. 
The GaussFit model has the ability to allow the 
offsets between data sets to be a free parame- 
ter. This is important because the radial veloc- 
ities from OAO and AAT are not absolute ra- 
dial velocities, but rather are measured relative 
to the iodine-free stellar template. Each data 
set thus has an arbitrary zero-point offset which 
must be accounted for in the orbit-fitting proce- 
dure (jWittenmver et al. I I2OO9I) . These offsets are 
3.8±3.2ms-i (OAO) and -0.1±3.7 ms'^ (AAT). 

For the final orbit fitting, we added radial- 
velocity jitter in quadrature to the velocity uncer- 



tainty of each observation. This jitter arises from a 
number of phenomena intrinsic to the star, such as 
granulation, o scilla tions, and activity (iSaar et al. 
[19981: IWrightli2005h . As HD4732 is an evolved 



star, scaling relations ( Kieldsen fc Beddind 19951 
20111 ) are not applicable. We thus determine the 



appropriate level of jitter empirically, by perform- 
ing two-Keplerian fits with varying levels of jit- 
ter. On some dates, the AAT recorded multiple 
consecutive exposures of HD 4732 (Table 13]). This 
provides one way to estimate the jitter: by noting 
the velocity spread of the multiple exposures. On 
the four dates which had multiple exposures, the 
mean velocity spread is 2.9m s~^. This is, how- 
ever, a lower limit as the timescales of many jit- 
ter sources are on the order of hours (granulation) 
or months (rotation). We also try a range of jit- 
ters from 4-10 ms^^, and note the results of the 
fits in Table 21 The estimation of radial- velocity 
jitter is quite uncertain, and this is reflected in 
Table HI where we see that the amount of added 
jitter has little effect on the quality of the result- 
ing fit. We adopt a jitter estimate of 6m s~^, as 
this value brings the reduced of the fit close to 
unity. The resulting 2-Keplerian planetary param- 
eters are given in Table [S] The uncertainties on 
the fitted parameters are derived in the usual way, 
from the covariance matrix of the least-squares fit. 
The phase gap for planet b has the greatest im- 
pact on the uncertainty in Tq for that planet. As 
a result, the uncertainty (Jt=18.4 days remains 
large even though many orbital cycles of planet b 
have been observed. The two-planet fit has a re- 
duced of 1.05, and the RMS scatter about the 
fit is 7.09ms~-^. Using a stellar mass of l-74lg'2o 
20081) . the planets have mini- 



M^I Takeda et al. 



mum masses m sin i of 2.4±0.3 (inner planet) and 
2.4±0.4 Mjup(outer planet). 

5. Line Shape Analysis 

In order to investigate other possible causes of 
observed radial-velocity variations such as pulsa- 
tion and rotational modulation rather than orbital 
motion, spectral line shape analysis was performed 
with use of high resolution stellar template spec- 
tr a. Details of the analy sis method are de scribed 
in lSato et afl (|2007t ) and lSato etall (|2002l) . 



At first, two stellar templates were extracted 
from five star-|-l2 spectra at phases with differ- 
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ent radial-velocity level, ~ 40 m s ^ and 



-20 



m s , of the observed radial velocities based on 



the method bv ISato et al. ( 2002 ). Cross correla- 
tion profiles of the two templates were calculated 
for about 90 spectral chunks (4-5A width each) in 
which severely blended lines or broad lines were 
not included, and then three bisector quantities 
were derived for the cross correlation profile for 
each chunk: the velocity span (BVS), which is 
the velocity difference between two flux levels of 
the bisector; the velocity curvature (BVC), which 
is the difference of the velocity span of the upper 
half and lower half of the bisector; and the velocity 
displacement (BVD), which is the average of the 
bisector at three different flux levels. Here we used 
flux levels of 25%, 50%, and 75% of the cross cor- 
relation profile to calculate the above three bisec- 
tor quantities, and obtained BVS=3.9±4.7 m s""'^, 
BVC=2.4±3.3 m s'^, and BVD= --57±11 m s^^ 
for HD 4732 (Figure O. Both of the BVS and the 
BVC are identical to zero and the average BVD 
agrees with the velocity difference between the two 
templates. Then the cross correlation profiles can 
be considered to be symmetric and thus the ob- 
served radial- velocity variations are best explained 
by parallel shifts of the spectral lines rather than 
distortion of them, which is consistent with the 
orbital-motion hypothesis. 

6. Dynamical Stability Analysis 

We here present dynamical studies of the HD 
4732 system. Since the eccentricities of the planets 
are relatively small and it is thus considered that 
they probably have not experienced close encoun- 
ters, we here assume that the both planets share 
the same orbital plane and are prograde. For the 
numerical integrations, we use a 4-th order Her- 
mite scheme. Figure |6] shows the one million year 
evolution of the best-fit two-Keplerian model de- 
rived in section HI for i = 90° orbits (i.e. planet 
masses at minimum). The system has two secular 
eigenfrequencies related with eccentricity, ll"/yr 
and 28"/yr. Accordingly, the planetary pericen- 
ter, a(l — e), and apocenter, a(l + e), distances 
oscillate with the secular timescale of 76000 yr. 

For lower inclinations, since the mass increases 
proportional to f/sini, mutual perturbations be- 
come stronger. The stability map in the (a, e) 
plane of the outer planet for different values of 



i is shown in Fig. [T] Here we plot the stabil- 
ity index D = \ (n2) — (ni)| (in °/yr) following 



when logj^Q D < —3 (blue to navy), while the sys- 
tem is chaotic when logj^Q D > —1 (yellow to red) 



Couetdic et al. I ( 2010l ). which studied stability of 
HD20220 6 system using Laskar's frequenc y map 
analysis (jLaskar I [l990l: iLaskar et al. l(200ll) . Fol- 
lowing the analysis, we calculate an average of 
mean motion of the outer planet in 1000 Kepler 
periods ((ni)) and subtract it from an average 
of mean motion of the same planet obtained in 
the next consecutive 1000 Kepler periods ((^2)). 
Through comparison with the five times longer 
simulations, we find that the system is regular 
D < 

D > 

As inclination decreases, the system loses stabil- 
ity. There is no stable area around la of the best 
fitted coplanar orbits for j < 5°. Thus, in the 
coplanar prograde configuration, we can roughly 
set the upper limit on the mass of both of the two 
planets to be about 28Mjup (i > 5°), which falls 
in the substellar mass regime. 

In top left panel of Fig. [3 one can find effects 
of mean-motion resonances as green vertical inci- 
sions. In the case of this system, based on the best- 
fit values of argument of pericenter, orbits closer 
to the mean-motion resonances are more unstable. 
The best-fit solution derived by the radial- velocity 
data (white cross) sits between the 7:1 and 8:1 
mean-motion resonances. From the orbital inte- 
gration shown in Fig. [6l we have also confirmed 
that their resonant angles are not librating. We 
thus conclude that the outer planet is not in the 
mean-motion resonances of the inner planet. 

7. Summary 

We report the detection of a double planet sys- 
tem around the evolved intermediate-mass star 
HD 4732 (KO IV, M = 1.7 M©) by precise radial- 
velocity measurements at OAO and A AO. The 
system is composed of two giant planets with min- 
imum masses of 7712 sini — 2A Mj, orbital period 
of 360.2 d and 2732 d, and eccentricity of 0.13 and 
0.23, respectively. The joint observations of OAO 
and AAO allowed us to increase the phase cov- 
erage of the nearly 1-yr periodicity of the inner 
planet and detect long-period signal of the second 
outer planet. 

The configuration of the system is similar to 
those of other multiple-planet systems currently 
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known: two giant planets in a combination of 
intermediate- and long-period orbits with rela- 
tively low eccentricities. HD 4732 c is the outer 
most planetary candidate (mpSini < 13 Mj) ever 
discovered around intermediate-mass stars except 
for those discovered by direct imaging The pe- 
riod ratio of the two planets is close to 7:1 or 8:1, 
but dynamical analysis showed that the system 
is not in the mean-motion resonance and orbits 
close to mean-motion resonances are unstable in 
this case. 

One million year orbital evolution of the best- fit 
two-Keplerian model in a coplanar prograde con- 
figuration showed no significant changes in eccen- 
tricities, but showed their slight oscillations with 
secular timescale of 76000 yr. Based on the dy- 
namical stability analysis, we found that the sys- 
tem is dynamically unstable in the case of « < 5° 
for a coplanar prograde configuration. Thus, we 
can set the upper limit on the mass of both of 
the two planets to be about 28 Afj {i > 5°) in 
this case, which falls well into the substellar mass 
regime. 
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Oph c (msini = 27 Afj, a=6.1 AU: ISato et al.ll2012fl 
is the outer most brown-dwarf companion ever discovered 
around intermediate-mass stars by radial velocity measure- 
ments. 
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Fig. 1. — HR diagram for HD 4732. Pairs of evolu- 
tionary tracks from Lejeune and Schaerer (2001) 
for stars with Z ~ 0.02 (solar metallicity; solid 
lines) and Z — 0.008 (dashed lines) of masses be- 
tween 1 and 3 Mq are also shown. 
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Fig. 2. — Spectra in the region of Ca H lines. Stars 
with similar spectral type to HD 4732 in our sam- 
ple are also shown. A vertical offset of about 0.3 
is added to each spectrum. 



Fig. 3.— Two-planet fit for HD 4732 radial veloc- 
ities. The planets have periods of 360 and 2732 
days, and the RMS about this fit is 7.09m s~^. 
OAO data are shown as filled black circles, and 
AAT data are filled cyan circles. 
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Fig. 4.— Left: Phase plot for HD 4732b (P = 
360 days), after removing the signal of the outer 
planet. Two cycles are shown for clarity. Right: 
Radial- velocity observations for HD 4732c (P = 
2732 days), after removing the signal of the inner 
planet. The symbols have the same meaning as in 
Figure m 
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Table 1 

Stellar Parameters for HD 4732 



Parameter 


Value 


Spectral Type 


KO IV" 


n (mas) 


17.70±0.99 


V 


5.90 


B-V 


0.944 


Av 


0.10 


Mv 


2.04 


B.C. 


-0.27 


%s (K) 


4959±25'' 


logs (cgs) 


3.16±0.08'' 


Vt (km s~^) 


1.12±0.07'' 


[Fe/H] (dex) 


O.OlzbO.04'' 


L (Lq) 


15.5 


R (Re) 


5.4 (5.0-5.8)'= 


M (Mq) 


1.74 (1.60-1.94)'= 


V sinirot (km s~^) 


1.45 


iTHip (mag) 


0.006 



''The star is listed in the Hipparcos cat- 
alogue as a KO giant. But judged from 
the position of the star on the HR dia- 
gram (Figure [T|, the star should be bet- 
ter classified as a less evolved subgiant. 

''The uncertainties of these values are 
internal statistical errors (for a given 
data set of Fe l and Fe ll line equivalent 
widths) evaluated by the procedure de- 
s crib ed in subsection 5.2 of lTakeda et al.l 
12002). 

"^The values in the parenthesis corre- 
spond to the range of the values assuming 
the realistic uncertainties in A log L cor- 
responding to parallax errors in the Hip- 
parcos catalog, AlogTcff of ±0.01 dex 
(~ ±100 K), and A[Fe/H] of ±0.1 dex. 
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Table 2 

OAO Radial Velocities for HD 4732 



JD-2400000 


Velocity (ms ) 


Uncertainty (ms ) 




-55.9 


4.7 


OOOO - U J- ^ t7 1 


-11.6 


3.7 




-57.7 


4.3 




-69.4 


6.4 




-62.8 


4.1 


^SfiQS QS447 


-30.2 


3.8 


^S7?S Q'nfiQT 


-4.8 


8.1 




-14.4 


4.3 


^'\7AA Q4fin7 


-3.3 


4.1 


Fi'^Qfi^ 2fi7R1 

t-j i_> c/ . ^ f cjx 


-64.7 


3.8 


Fi4mS 1 1Q7^ 


-37.6 


3.8 


iJ^U'i J . U ^ij O J. 


-17.8 


4.4 




-8.3 


3.3 




3.0 


6.2 




-47.9 


3.8 


^4^"^^ 2S2S8 


-38.5 


4.4 




-22.6 


4. 1 


^441 P, 09081 


-18.8 


4.8 


Fi44fin Q'^7Q'^ 


24.1 


3.0 




33.7 


9.5 


54496.90370 


46.0 


4.9 


54703.25119 


-38.0 


3.6 


54756.10553 


-7.5 


3.5 


54796.00048 


23.5 


3.8 


54817.94637 


31.3 


3.5 


55101.16674 


-5.9 


4.0 


55108.20638 


2.7 


4.6 


55133.10162 


11.9 


3.8 


55137.09148 


2.9 


4.1 


55161.03266 


32.6 


3.7 


55162.00982 


34.0 


2.8 


55183.92889 


42.0 


3.6 


55203.89392 


46.3 


3.5 


55233.90048 


68.5 


5.8 


55399.28930 


-13.2 


3.7 


55438.17586 


-22.0 


4.4 


55444.21320 


-29.9 


3.5 


55468.16146 


0.7 


5.0 


55480.13068 


13.1 


4.3 


55503.02690 


28.3 


3.9 


55510.06421 


22.8 


3.7 


55524.97523 


39.5 


3.4 


55546.92022 


41.9 


3.6 


55579.91720 


53.6 


6.0 


55787.29129 


-14.6 


4.8 


55877.01304 


6.5 


3.5 


55919.93294 


29.2 


4.0 


55937.88506 


41.8 


3.4 
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Table 3 

AAT Radial Velocities for HD 4732 



JD-2400000 


Velocity (ms ) 


Uncertainty (ms ) 


55455.16094 


-19.7 


5.6 


55455.20587 


-20.5 


5.9 


55456.06022 


-24.3 


5.5 


55494.91653 


21.9 


4.0 


55579.92330 


63.1 


4.8 


55600.90903 


67.7 


5.9 


55600.91897 


75.7 


5.1 


55600.93330 


71.8 


5.3 


55707.30755 


-3.1 


4.3 


55707.31526 


-4.4 


4.4 


55760.15572 


-36.3 


4.0 


55783.13308 


-37.3 


6.9 


55783.14079 


-35.9 


5.8 


55841.99044 


-12.3 


4.3 


55878.95827 


1.0 


4.4 


55880.01556 


10.3 


4.5 


56060.29762 


-31.2 


5.6 


56089.30186 


-39.6 


4.2 


56091.27740 


-46.8 


4.2 



Table 4 
Jitter Estimates 



Added jitter (ms ^) 


xl 


RMS (ms-^) 


2.9 


2.12 


7.10 


4.0 


1.66 


7.09 


5.0 


1.32 


7.09 


6.0 


1.05 


7.09 


7.0 


0.85 


7.08 


8.0 


0.70 


7.08 


9.0 


0.59 


7.08 


10.0 


0.49 


7.08 



Table 5 
Keplerian Orbital Solutions 



Planet 


Period 


To 


e 




K 


m sin i 


a 




(days) 


(JD-2400000) 




(degrees) 


(ms-i) 


(Mj„p) 


(AU) 


HD 4732 b 


360.2±1.4 


54967±18 


0.13±0.06 


85±16 


47.3±3.5 


2.37±0.34 


1.19±0.05 


HD 4732 c 


2732±81 


56093±103 


0.23±0.07 


118±15 


24.4±2.2 


2.37±0.38 


4.60±0.23 
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Fig. 5. — Bisector quantities of cross-correlation 
profiles between the templates of HD 4732 at 
phases with ~ 40 m s~^ and with ~ — 20 m s""*^ 
of the observed radial velocities: bisector veloc- 
ity span (BVS, circles), bisector velocity curvature 
(BVC, triangles), and bisector velocity displace- 
ment (BVD, squares). Offsets of 300 m s^^ and 
—300 m s^^ are added to BVS and BVD, respec- 
tively. 
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Fig. 6. — Evolution of semimajor axis a, pericen- 
ter a(l — e), and apocenter a(l + e) distances for 
the two planets of HD 4732 system for the best-fit 
two-Keplerian model to the radial-velocity data. 
Red lines are for semimajor axis and blue ones 
are for pericenter (lower line for each planet) and 
apocenter (upper line for each planet) distance. 
We here assume i = 90° and prograde coplanar 
orbits. 
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Fig. 7. — Stability map of the outer planet (HD 
4732 c) for a coplanar prograde configuration with 
different values of inclination i. Orbital elements 
other than semimajor axis and eccentricity of the 
outer planet are taken from the best fit values. 
The color scale shows the level of mean motion 
diffusion, logj^Q D. Blue to navy orbits are stable, 
while yellow to red orbits are chaotic. The step 
size of eccentricity is Ae=0.005. The step sizes of 
semimajor axis are Aa=0.005 AU (top-left panel), 
Aa^O.Ol AU (top-right panel), and Aa=0.05 AU 
(bottom panels). White crosses represent the best 
fitted a and e with their 1 a errors. Black dotted 
line is the orbit crossing line. 
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